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ABSTRACT 



We present the first survey of resolved stellar populations in the remote outer halo of our 
nearest giant elliptical (gE), Centaurus A (D= 3.8 Mpc). Using the VIMOS/VLT optical cam- 
era, we obtained deep photometry for four fields along the major and minor axes at projected 
elliptical radii of ~ 30 - 85 kpc (corresponding to ~ 5 - 14R e ff). We use resolved star counts 
to map the spatial and colour distribution of red giant branch (RGB) stars down to ~ 2 mag- 
nitudes below the RGB tip. We detect an extended halo out to the furthermost elliptical radius 
probed (~ 85 kpc or ~ 147? e ff), demonstrating the vast extent of this system. We detect lo- 
calised substructure in these parts, visible in both (old) RGB and (intermediate- age) luminous 
asymptotic giant branch stars, and there is some evidence that the outer halo becomes more 
elliptical and has a shallower surface brightness profile. We derive photometric metallicity 
distribution functions for halo RGB stars and find relatively high median metallicity values 

(<[Fe/H]> me j 0.9 to -1.0 dex) that change very little with radius over the extent of our 

survey. Radial metallicity gradients are measured to be « -0.002 - 0.004 dex/kpc and the 
fraction of metal-poor stars (defined as [Fe/H]< -1.0) is « 40 - 50% at all radii. We discuss 
these findings in the context of galaxy formation models for the buildup of gE haloes. 

Key words: galaxies: photometry - galaxies: stellar content - galaxies: individual: Cen A - 
galaxies: evolution 



1 INTRODUCTION 

Giant elliptical (gE) galaxies contain a large fraction of the stellar 
mass in the Universe and it is crucial to understand how they form 
and evolve. Important clues come from analyses of their structures, 
stellar populations and kinematics and a variety of models have 
been invoked in order to reproduce and explain their observed prop- 
erties. The classical monolithic collapse models have stars form- 
ing during an early dissipative phas e where the gas rapidly flows 
into the central galaxy regions (e.g. lLarsonlll974l: ICarlberdfl984l : 
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lArimoto & Yoshiilll987h . On the other hand, in the hierarchical 
scenario of a ACDM cosmology, either a few major mergers or 
several minor mergers and accretion event s are believed to build- 
up the haloes of the largest galaxies (e.g. iKauffmann et alj[l993l : 
iDe Lucia et al.ll2006[). Recent work has favoured a two-stage sce- 
nario (e.g. Naab et al.ll2007L 120091 : iKavirai et al]l2009l : bser et all 
I2O10L 2012) whereby early in- situ star formation builds up the cen- 
tral parts of gEs at z > 2, while later low-mass dry accretions feed 
the growth at large radii. Such models provide a natural explana- 
tion for the observation that some massive ellipticals at high z are 
more compact than their low redshift counterpa rts (e.g.. lDaddi et al.l 
l2005l : lTruiillo et al.ll2006l : lBuitrago et aDl2008h . 

Observations confirm the existence of a dominant popula - 
tion of old stars in the centers of gEs (e.g. Thom as et al .2005), 
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but few constraints exist to date on the properties of stars in re- 
gions beyond the effective radius, R eS . The reason for this is that 
most studies of gEs have been based on integrated light analy- 
ses, and the faintness of the outer regions has hindered quantitative 
study of these parts. Some work has demonstrated the existence of 
faint tidal debris in the outskirts of ellipticals (e.g. iTal et al .1 12009; 
Ijanowiecki et al.1 EoiOh while a few recent studies have started 
to probe kinematics and chemical composition out to increa s- 
ingly large radii (e.g . Jweii mans et al.ll2009l:ISpolaor et al.ll2010allbl ; 
ICoccato et al.l2010ai lbF lGreene et alfcOllLa Barbera et al.l2012l) . 
The general picture emerging is one of older ages and lower metal- 
licities in the outskirts of gEs, qualitatively consistent with the pre- 
dictions of the two-phase model. However, it should be noted that 
these studies typically have not probed beyond a few effective radii. 

To probe gEs to even more extreme rad ii, techniques other 
than integrated light analysis are required. & van Dok kum 
(201 J) stack SDSS images of more than 42000 massive red galax- 
ies at z ~ 0.35 in order to explore the mean properties to radii 
of ~ 400 kpc. They find the halo colour is bluer than the central 
regions and remains constant beyon d ~ 3R^ . The y also find a de- 
viation from a simple Sersic profile (Sersic 1968)) at radii > 8/? e ff» 
in the sense that there is excess light, and an increasing elliptic- 
ity. While powerful, the stacking technique has the disadvantage 
that information regarding individual galaxies is lost and the in- 
terpretation can be complicated if the stacked sample is not truly 
homogeneous. 

An alternative method is to use resolved stellar populations to 
probe the outskirts of gEs. This technique can only be applied to 
a small number of very nearby systems at present, but has the ad- 
vantage that age and metallicity can be directly constrained from 
colour-magnitude diagram (CMD) morphology, removing several 
of the uncertainties inherent in integrated light analyses. Wide-area 
resolved star counts also allow the structure and extent of individual 
galaxies to be probed to extremely low effective surface brightness 
levels, competitive with those attained through the stacking analy- 
sis of many thousands of systems. This has already been demon- 
strated in a number of recent studies of spiral and dwarf galax- 



ies (e . g.lFerguson et al.l2002l;llbata et alJ :2007 ; 



2009; Barker et al. 2009, 120121; iBernard et al. 



McCon nachie et al] 



120121) but thus far 



there has been no wide-field study of the resolved stellar popula- 
tions in a gE. 

Located at a distance of 3.8 Mpc (lHarris et"al1l2010ah and 
sitting at the centre of the homonymous group of galaxies, Cen- 
taurus A (Cen A, NGC 5128) is the best target for wide-field re- 
solved stellar populations analysis of a gE from the ground. This 
peculiar system hosts an active nucleus and shows striking evi- 
dence for having experienced a recent gas-rich merger event. In 
particular, the inner regions (< 30 kpc, or within the 25th 5-band 
isophote, R25) contain dust lanes, a dense warped gas disk and dif- 
fuse stellar shells and arcs (e.g. lHavneseUd] [l983; MalinetalJ 
1983; |Peng etaflbOOl see also Fig. Q). The strong radio emis- 
sion from Cen A partly overlaps with these optical features, but 
extends very much furthe r out, c overing almost 8 degrees in dec- 
lination (for a review, see Israel 1998). None t heless, with a mass 
of ~ .5 - 1 x 10 12 M Q dWoodlev et al.ll2007l ; IWoodlev & Gomezl 
2010c) and a luminosity of My ~ -21.5, Cen A is typical in a 
global sense of the field ellipticals seen at low redshift. 

Fifteen years ago, the old stellar populations of Cen A were 
the first to be resolved in a galaxy outside the Local Grou p by using 
WFP C2 onboard the Hubble Space Telescope (HST) dSoria et all 
1996). Since then, numerous space- and ground - based studies 
(lHarris et al] Il999l; lHarris & Harris] l2000l b002l iMarleau et al.l 



120001 ; iReikuba et alfeOOll. 120031 . l200l l201ll : ICrockett et al]l2012h 
have surveyed Cen A in both optical and near-IR (NIR) bands 
out to radii of ~ 38 kpc in projection. At the distance of Cen A, 
the physical scale i s 1 arcmin ~ 1 kpc and, given that R eS ~ 
330 arcsec~ 6.1 kpc (I van den B ergh 1976), this radial extent corre- 
sponds to ~ 6.2/?efF- Amongst the results, the predominantly old red 
giant branch (RGB) population is found to be moderately metal- 
rich ([M/H]~ -0.65) with a broad metallicity spread and no sig- 
nificant metallicity gradient. In the deepest HST field studied to 
date (~ 38 kpc), an intermediate- age population is estimated to ac- 
count for ~ 20 - 30% of the p opulation, and has a derived age of 
~ 2-4 Gyr feeikuba et al.l201ll) . A much younger (~ 10 Myr), and 
more localized, trace population is found within ~ 8 kpc, aligned 
with the radio jet of Cen A and possibly trig gered by bow shock s 
propagating through the interstellar medium (ICrockett et al.ll2012h . 
Although probing a significant range in radius, these studies are all 
based on small field-of-view (FOV, a few arcmin 2 ) imagery, leav- 
ing open questions about how representative they are of the global 
properties of Cen As stellar halo. 

There is also good reason to believe that the halo of Cen A may 
extend very far out, well beyond the region where previous resolved 
stellar populations studies have probed. Indeed, globular clusters 
(GCs) and planet ary nebulae have b ee n discovered out to projected 
radii of ~ 85 kpc jPeng et dfeOCMafllWoodlev et al.ll2007lJ2010bl ; 



lHarris et alfcoil see also Fig. HI. lHarris et al.l d2007ah discuss the 
halo of another nearby gE, NGC3379, where the outer stellar halo 
is found to be predominantly metal-poor at radii > 10^ e ff • Such a 
transition from a metal-rich to a metal-poor halo has not yet been 
observed in Cen A suggesting that, if it is a universal feature of gEs, 
it may lurk at larger radii. 

With this motivation in mind, we have undertaken the first 
wide-field survey of the resolved stellar populations in the remote 
outer halo of Cen A. This work is part of a larger programme we 
are conducting to explore the low surface brightness outer regions 
of galaxies within 5 Mpc using wide-field imagers on 8-m class 
telescopes (iBarker et al.l2009Ll2012l;lBernard et al.b012b . In g^we 
present our observations and photometric analysis, and in |3] we 
present the resulting CMDs. In §|4] we investigate the structure of 
Cen As outer halo while in |5]we derive metallicity distribution 
functions (MDFs) and constrain the radial gradients. In |6]we dis- 
cuss our results, and the summary is presented in |7] 



2 THE DATA 
2.1 Observations 

Observations were obtained in imaging mode with the VIMOS in- 
strument at the ESO Very Large Telescope (VLT) on Cerro Paranal, 
Chile under programme 074 B -0741 (Cycle 74, PI: Ferguson). VI- 
MOS (Le Fevre etaD l2003h is a visible wide field imager and 
multi-object spectrograph mounted on the Nasmyth B platform of 
the 8.2m Unit Telescope 3 (Melipal). The instrument is made up 
of four identical arms, each with a field of view of 7 x 8 arcmin 2 , 
separated by ~ 2 arcmin. The pixel scale is 0.205 arcsec. 

We imaged two fields along each of the major and minor axes, 
centered at projected galactocentric elliptical radii of R ~ 45 and 
70 kpc; these fields all lie well beyond the the known shell system 
of Cen A which is generally confine d to < R25 (lHavnes et al.ll 19831 ; 
iMalin et al.lfl983l ; IPeng et al.ll2002h . The radii, which refer to the 
semi-major axis of the elliptical isophote which intersects a given 
position, are calculated assuming an axis ratio of b/a= 0.77 and a 
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Figure 1. Projected position on the sky in standard coordinates (with respect to the center of Cen A) of the four VIMOS observed fields (black rectangles). 
Left panel. We show a photographic plate image of the central regions of Cen A. The solid red ellipse represents the 25th 5-band isophote as listed in the 
RC3 (R25 = 25.7' ~ 27.8 kpc, b/a= 0.77 and PA= 35°), while red dashed ellipses are drawn at projected radii of 45 kpc (corresponding to ~ 7.5/? e ff)> 
65 kpc (~ lO^eff) and 85 kpc (~ 14R e ff)- The black solid lines represent the major and minor axes. The blue sq uare shows the location of the outermost 
ACS/HST pointing (projected radius of ~ 38 kpc) in the Cen A halo to d ate, analysed b y Reikuba et al. ( 20051 12OI1I) . Right panel We show the positions of 
confirmed globular clusters (red circles) from the IPene etalJ J2004al) and lWoodlev et alj (2010b) catalogs, and confirmed planetary nebulae (blue asterisks) 
from lPeng et al.l (|2004b). The latter exten d out to projected rad ii of 85 kpc. The two magenta stars indicate the positions of two early-type dwarf companions 
of Cen A (KK197 and KKs55, studied in lCrnoievic etalkoiCh . 



Table 1. Observing log for each of the four pointings. 



Field 


O72000 


£72000 
(° ' ") 


T? a 

(arcmin/kpc) 


Kit 

(arcmin/kpc) 


F 


(sec) 


N exp 


Seeing med 
(") 


Cen A-maj 1 


13 27 39.2 


-42 26 35.2 


42.3/45.7 


42.3/45.7 


/ 


6000 


15 


0.73 




13 27 40.7 


-42 26 30.1 






V 


10000 


8 


0.76 


Cen A-maj 2 


13 28 55.1 


-42 09 57.6 


63.8/68.9 


63.8/68.9 


I 


5200 


13 


0.68 




13 28 55.0 


-42 10 03.7 






V 


10000 


8 


0.58 


Cen A-minl 


13 27 49.8 


-43 20 13.9 


32.2/34.8 


41.8/45.1 


I 


8000 


20 


0.57 




13 27 49.8 


-43 20 17.2 






V 


8750 


7 


0.66 


Cen A-min2 


13 29 04.9 


-43 31 29.3 


49.8/53.8 


64.7/69.9 


I 


6800 


17 


0.68 




13 29 05.0 


-43 31 33.6 






V 


10000 


8 


0.60 



a Projected linear distance from Cen As center (072000 = 13 /z 25 m 27.6 s , £72000 = -43°01'08.8"; taken from the NED 
database). 

^ Projected ellipti cal radius from Cen As center, computed assuming b/a= 0.77 and PA= 35° (values from the RC3; 
Corwin et al. 1994). 



position angle of 35° JCorwin et al.lll994h . Throughout this paper, 
whenever we refer to a "radius", we mean a projected elliptical ra- 
dius computed in this way. The footprint of each pointing is drawn 
and labelled in Fig. [T] where some fiducial radii are also marked. 
This figure also shows the location and size of the deepest extant 
HST/ACS pointing which lies at a projected radius of ~ 38 kpc. 
Each of our VIMOS fields was observed with 13 to 20 individ- 
ual exposures of ~ 400 sec in the /-band, and 7 to 8 exposures 
of ~ 1250 sec in the V-band. Individual exposures were dithered 
by a few tens of arcsecs. The seeing was excellent during most of 
the observing runs, ranging from 0.45-0.98" with median values of 



~ 0.67" for /-band images and ~ 0.65" for V-band images. Full 
details of the observations are reported in Tab. [T] 

The sensitivity of each individual VIMOS detector is lower at 
the edges, with additional distortion and vignetting in some cases. 
Stars which fall in these regions are thus not recovered well and 
we exclude these regions from our analysis. The effective area that 
was surveyed by our observations (i.e., after eliminating inter-chip 
gaps, masking bad edges, and considering the overlap between the 
minor axis fields) is ~ 0.22 deg 2 , corresponding to ~ 30 x 30 kpc 2 . 

The data reduction procedure followed a modified version of 
the pipeline developed for processing Wide Field Cam era data from 
the Isaac Newton Telescope (for further details see Ilrwinlll985l 
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19971: llrwin & Lewisll200ll : llrwin et all Eooi barker et all E0T2I : 
Bernard et al J 12012 ). The first stage includes bias and overscan- 
correction, plus trimming of each image to the useful active detec- 
tor area. Master flats were then created by stacking a well-exposed 
dithered set of 21 V-band and 9 /-band twilight sky observations. 
The flat- fielding also corrects for internal gain differences between 
the detectors. The dark sky /-band images were combined to form a 
master fringe frame. Typical fringe pattern amplitudes were around 
4% of the sky with a range of spatial scales. Individual /-band sci- 
ence images were corrected for fringing using a scaled version of 
the master fringe frame. 

Before stacking the science images, object catalogues were 
generated for each image to refine the astrometric calibration 
and to assess the data quali ty. A Zenithal polynomial projection 
dGreisen & C alabretta Eo02h ' was found to provide an accurate 
World Coordinate System (WCS) of the telescope plus imaging 
camera. A third-order polynomial included all the significant radial 
field distortions and a further six-parameter linear model per detec- 
tor was used to define the remaining astrometric transformations. 
The Two Micron All Sky Survey (2MASS) point-source catalogue 
dCutri et al.ll2003l) was used for the astrometric reference system. 

Common background regions in the overlap area of individ- 
ual exposures were used to correct for sky variations during the 
exposure sequence and the final stack included seeing weight- 
ing, confidence (i.e., variance) map weighting and clipping of cos- 
mic rays. Fi rst-pass catalogue s using standard aperture photometry 
techniques (llrwin et al .1120041) were generated and used to update 
the astrometric solution for the stacked images and to provide the 
aperture corrections to total flux required for the final photometric 
calibration. 

I n addition to th e science targets, a series of Landolt standard 
fields (Landolt 1992) were observed throughout the programme. 
These were analysed in the same way as the science image and 
provide the instrumental to photometric system (Johnson-Cousins) 
calibration. The variation in derived magnitude zero-points seen in 
the standard field observations was at the level of +2% and we take 
this to be indicative of the final accuracy of our photometric cali- 
bration. 



2.2 Photometry 

Although aperture photometry was performed as part of the data 
reduction process, we chose to further perform point spread func- 
tion (PSF)-fitting photometry on the stacked images so as to more 
accurately recover faint stars in crowded regions, obtain additional 
morphological information on detected sources and better assess 
photometric errors and incompleteness effects. For this we used 
the photometric suit e of programs DAOPHOT and ALLFRAME 
dStetsonll 19811 1994b . For each chip of the stacked images, the PSF 
was constructed by using at least 30 bright, non saturated stars, 
evenly spread across the image. We derived a list of objects with 
a 3cr detection threshold in each filter, and computed the coordi- 
nate transformations between the different filter s with the DAO- 
MATCH/DAOMASTER packages dStetsorJ[l993h . The list of stars 
detected in at least one of the two filters for each chip was used as 
the input list for the simultaneous photometry on the stacked im- 
ages using ALLFRAME. In the final stellar catalogue, we retain 
sources that were detected in both filters with photometric errors 
smaller than 0.3 mag. Furthermore, we only retain sources that are 
highly likely to be stellar by requiring absolute values of the sharp- 
ness parameter < 2. At magnitudes fainter than / ~ 23 , we also 



Table 2. /o-band 50% completeness levels at different colours and elliptical 
radii. 



(V-I)o 


0-1 


1-2 


2-3 


R< 45 kpc 


Cen A-maj 


25.2 


24.8 


24.2 


Cen A-min 


25.7 


25.1 


24.3 


R> 45 kpc 


Cen A-maj 


25.4 


25.0 


24.4 


Cen A-min 


25.6 


25.2 


24.4 



require the^f parameter to be smaller than 1.4-1.6, depending on 
chip. 

We calibrate the magnitudes obtained with our PSF photom- 
etry to the Johnson-Cousins system /- and V-bands by selecting 
bright stars in common with the aperture photometry catalogue, 
and then deriving a constant offset between the two. The total num- 
ber of stars remaining after quality cuts is 104325, of which 37117 
are in the inner major axis field (Cen A-maj 1), 22170 are in the 
outer major axis field (Cen A-maj2), 26922 are in the inner minor 
axis field (Cen A-minl), and 18046 are in the outer minor axis field 
(Cen A-min2). There is partial overlap between one chip in each of 
the Cen A-minl and Cen A-min2 pointings. We chose to match the 
common sources and retain the photometry from the Cen A-minl 
field, since this is the one with the longest exposure time and best 
seeing in the /-band, hence smaller photometric errors. 

2.3 Artificial Star Tests 

We have performed extensive artificial star tests to quantify the ob- 
servational uncertainties in our data. We have simulated fake stars 
for each individual chip in each field with an even spatial distribu- 
tion and with colours and magnitudes covering the entire colour- 
magnitude space spanned by the observed stars. The total number 
of simulated stars per chip is ~ 250000, divided into ~ 80 sub- 
sets so as to not substantially increase the stellar crowding on the 
images. The artificial stars were added to each chip by using the 
measured PSF. The photometry was then performed in the same 
way as for the real data, with the same quality cuts applied. We de- 
rive photometric errors as the difference between the input and the 
recovered magnitudes of the fake stars, and assess the incomplete- 
ness of our data by counting the number of stars recovered over the 
total number of injected stars. 

Given that our pointings cover a significant portion of Cen As 
outer halo, it is important to quantify the photometric uncertainties 
not only as a function of magnitude but also as a function of galac- 
tocentric radius. In Fig.|2]we show the completeness curves for the 
entire sample (i.e., without distinction between major and minor 
axes) within and beyond 45 kpc and for various colour ranges. As 
expected, the completeness at a fixed magnitude is somewhat lower 
at smaller radii because of increased crowding, and it decreases to- 
wards redder colours. The derived / -band completeness values for 
the individual major and minor axes are given in Tab. |2] The radial 
completeness at a fixed colour does not vary significantly along the 
minor axis, partly due to the fact that the Cen A-minl pointing had 
a longer /-band exposure time than other pointings. On the other 
hand, the radial dependence of the completeness along the major 
axis is mainly driven by the intrinsically higher stellar density at 
smaller radii. Along the major axis, the typical photometric errors 
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Figure 2. /o-band completeness curves for the four VIMOS fields com- 
bined (i.e., without distinction between minor and major axes). The differ- 
ent curves refer to distinct colour ranges and elliptical radii, as indicated. 
The dash-dotted red line denotes the 50% completeness level. 



in magnitude are of order ~ 0.1 mag at I = 22.9 for radii < 45 kpc, 
and at Iq = 23.4 for R > 45 kpc. Along the minor axis, the pho- 
tometric errors are ~ 0.1 mag at a I ~ 23.7 - 23.8 for all radii. 
This information will be incorporated into our spatially-resolved 
analysis of Cen A's stellar populations whenever necessary. 



- +0 .2^ stellar isochrones from the Dartmouth group (iDotter et al.l 
l20Q8h shifted to the distance of Cen A with a fixed old age (12 Gyr) 
and metallicity varying from [Fe/H] = -2.5 to -0.3. For our choice 
of Qf-enhancement, these values translate into a range of [M/H] 

from 2.4 to -0.15 dex using the prescriptions of ISalaris et al.l 

dl993h . This set of theoretical models has been shown to give a par- 
ticul arly good fit to old and intermediate- age star cluster fiducials 
(e.g. lGlatt et alj2008p. We also show the AGB phase of the Padova 
stellar isochrones (iGirardi et al .1120 1 Ol) . We choose isochrones with 
a metallicity of Z = 0.004 (corresponding to [M/ H]- -0.7) and 
ages of 2 and 4 Gyr (following the age estimates of iReikuba et al.l 
l201ll for fields at smaller radii). 

Fig. H] also shows the selection boxes used to define the loca- 
tions of the RGB and AGB candidates that we will subsequently 
analyse. For the RGB box, the faint limit is set by the ~ 50% com- 
pleteness limit for radii > 45 kpc while the bright limit is deter- 
mined by considering the location of the TRGB. The blue end of 
the RGB box is defined by the most metal-poor isochrone and ex- 
cludes the bulk of unresolved contaminants. For the AGB box, we 
start ~ 0. 1 mag above the TRGB in order to avoid RGB stars scat- 
tered up by photometric errors and allow for a luminosity width of 
1.5 magnitudes. The colours range from the bluest edge of the RGB 
to (V - I)o ~ 2.8, thus not encompassing the reddest stars in the 
CMD. Stars with colours redder than this might be dust-enshrouded 
carbon-rich AGB stars, however only a handful of these are likely 
to be present and their numbers will be dwarfed by the foreground 
contaminants which dominate at these colours. 



3.1 The Distance to Cen A 



3 COLOUR-MAGNITUDE DIAGRAMS 

The CMDs for each of the four VIMOS pointings are shown in 
Fig. [3] We have deredde ned the original m agnitudes using the 
Schlegel extinction maps dSchlegel et al.lll998h and adopting the 
ICardelli et al.l (|_1989) extinction law. In order to account for possi- 
ble spatial variations in the reddening (for example, due to the large 
FOV and the low Galactic latitude of Cen A), we have applied ex- 
tinction corrections to subfields of 20 arcsec on a side. On these 
scales, the variations are on the order of ~ 0.02 mag with respect 
to the mean extinction values of A/ ~ 0.22 and A v ~ 0.38. Internal 
extinction is expected to be negligible in the gas-deficient outer re- 
gions of Cen A and we do not take it into account in the subsequent 
analysis. 

The main feature visible in the CMDs of Fig. [3] is a promi- 
nent RGB, a feature indicative of an old population (> 1 Gyr) of 
evolved cool giants. The upper two magnitudes of the RGB can 
be seen in all cases. The width of the RGB is broader than the 
photometric errors derived from the artificial star tests implying an 
intrinsic spread in colour, likely due to metallicity. The absence of 
upper main sequence and blue or red supergiant stars excludes the 
presence of stars younger than ~ 1 Gyr. In the CMDs, there are 
also some stars brighter than the tip of the RGB (TRGB) which 
are candidate intermediate- age (~ 1 - 8 Gyr) luminous asymptotic- 
giant branch (AGB) stars. However, as discussed below, this region 
of the CMD is also heavily contaminated by foreground Galactic 
stars. Similarily, we will also show that objects in the colour range 
0.3 < (V - I)o < 1 are for the most part unresolved background 
galaxies. 

In the top panel of Fig. [4] we show the combined CMD for all 
stars in our photometric catalogue. We overlay ar-enhanced ([ar/Fe] 



The TRGB has a fixed /-band absolute magnitude of Mj_ 4.05 

for p redominantly metal-poor old populations (e.g., iRizzi et all 
l2007h and can therefore be used as a powerful distance indica- 



Harris et al 1 d2010ah reviews Cen As distance measurements 



tor. 

to date, as derived with different indicators. Based on their previ- 
ous HST studies, they report a mean Itrgb = 24.10 + 0.1 which, 
when combined with their assumed A/ = 0.22 + 0.02, yields 
Io,trgb = 23.88 + 0.1. We recompute the TRGB value s for each 
of ou r observed fields using a Sobel edge-detection filter dLee et al.l 
1 19931) and average them to find Iqjrgb = 23.81 + 0.35, in ex- 
cellent agreement with previous estimates. The significant uncer- 
tainty on our value comes from a combination of zero-point and 
extinction uncertainties, and our chosen binwidth (which depends 
on the photometric error at the TRGB value). The precise absolute 
magnitude of the T RGB has a colour/metallicity dependence (e.g. 
Mad ore et al] 12009), becoming fainte r for metal-rich p opulations 
([Fe/H]> -0.5). The value reported bv lHarris et aDfcOlOah is com- 
puted considering only the metal-poor side of the TRGB. Our data 
do not reach the reddest magnitudes of Cen A's RGB ((V - 1) > 3) 
because of incompleteness, so we too are sampling a region where 
the TRGB magnitude is mostly constant. To be conservative, we 
used stars with (V - 7) < 2.5. Finally, we exclude stars bluer than 
(V - I)o ~ 1 for the TRGB computation, as these are largely con- 
taminated by unresolved background galaxies. Throughout the pa- 
per we will use the TRGB value derived in this study. 



1 The choice of the ^-enhancement is justified by the mean value found for 
Cen A GCs, [a/Fe]~ +0.14 ( Woodlev et al. 20 njl). 
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Figure 3. Dereddened CMDs of each individual VIMOS field along the major and minor axes (labeled as in Fig. [I]. Photometric errorbars are reported, as 
computed from artificial star tests. The red line at Io = 23.81 aside each CMD indicates the average TRGB value derived from the four fields (see text for 
details). 



3.2 Background and Foreground Contaminants 

Even under excellent seeing conditions, a substantial number of 
high-redshift background galaxies will appear as unresolved point 
sources in ground-based data. Moreover, the low Galactic latitude 
of Cen A (b ~ 20°) means that a non-negligible number of Milky 
Way dwarf stars will be projected along our sight lines. Minimizing 
this contamination, and ultimately correcting for it, are necessary 
before we can properly interpret our data. 

The mini field was observed with the best combination of 
depth and seeing of all our fields. Fig. [3] shows that this CMD has 
two distinct almost parallel sequences, the particular definition of 
which is due to the smaller photometric errors. The blue sequence 
lies blueward of plausible RGB tracks (see top panel of Fig. [4]), 
but redward of where we would expect to see blue plume stars 
(V - I ~ 0). Similarly, red and blue supergiants at the distance 
of Cen A would be expected at slightly redder and bluer colours 
respective ly but at brighter m agnitudes than the sequence seen (see 
Fig. 9 of lBarkeretaLll2012h . As in previous work, we therefore 
identify t he blue sequence as resulting from unresolved b ackground 
galaxies (Barker et al. 2009, 2012; Bernard et al. 2012). The fact 
that the surface density of objects in this region of the CMD does 
not show any trend with galactocentric radius adds further weight 
to it being a population unassociated with Cen A. To further sup- 



port our interpretation, we additionally show in the middle panel 
of Fig. [4] the CMD for all the sources rejected as stars by quality 
cuts in our PSF photometry. These sources will largely be resolved 
galaxies. As can be seen, the dominant sequence here does indeed 
lie blueward of the main RGB sequence. We have used this fact to 
guide our choice of RGB selection box boundaries so as to limit the 
contamination by these non- stellar sources. While it is clearly im- 
possible to cleanly separate stars from unresolved galaxies through 
color-magnitude cuts alone, our conservative approach is designed 
to optimise the fidelity of our RGB catalogue, if not its complete- 
ness. 

To quantify the contamination from the Milky Way fore- 
ground, we have simulated the CMD of the Galactic pop ulation in 
this direction using the Besancon model (Robin et al. 2003). As in- 
puts, we provide the Galactic coordinates and approximate areas of 
each of the four observed VIMOS fields. We apply the photometric 
errors derived from our artificial star tests to the simulated stars, 
and account for incompleteness by randomly extracting a sub sam- 
ple reflecting the derived recovery level in each magnitude bin. The 
bottom panel of Fig. |H shows the simulated Galactic CMD for the 
four pointings combined. The densest sequences here are clearly 
recognisable in the observed CMDs and only minimally contam- 
inate our RGB selection box. However, the Galactic populations 
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severely contaminate the AGB selection box, where they contribute 
up to ~ 50% of the sources along the minor axis. Overall, the pre- 
dicted number of foreground stars accounts for ~ 10 - 15% of the 
total number of stars observed. 

Consideration of the magnitudes and colours of the back- 
ground and foreground contaminants therefore leads us to conclude 
that our analysis of Cen A's remote RGB population will not be 
strongly affected by their presence. In the subsequent sections, we 
will use two methods to estimate the actual contaminant source 
density in our RGB selection box. In the first instance, we will use 
knowledge of the surface density of objects rejected as stars (i.e. 
that are resolved) that fall within the RGB selection box to provide 
a lower limit on the contaminant level. We will also use the surface 
density of RGB stars in the 75 < R < 85 kpc radial bin along the 
minor axis as an alternative estimate of the contaminant level. This 
region has the lowest point source density of all those studied in 
our survey. Since the number of genuine Cen A RGB stars which 
reside in this area is unknown, it strictly provides an upper limit on 
the contamination. The situation is more difficult for the AGB se- 
lection box and we resort to using model predictions for the surface 
density of Galactic stars to assess the contaminant level here. 



4 HALO EXTENT AND STRUCTURE 

4.1 Colour-Magnitude Diagrams as a Function of Radius 

In Fig. [3] we show the CMDs for each of our individual point- 
ings; we now consider the CMDs as a function of radius in or- 
der to search for radial variations. Fig. |5]shows Hess diagrams for 
three radial bins along each of the major and minor axes with Dart- 
mouth stellar isochrones overlaid. Stellar density decreases with ra- 
dius along both axes, with the highest density measured anywhere 
occurring at R < 45 kpc along the major axis. A clear RGB se- 
quence can be seen at all distances along the major axis, testifying 
to the existence of a halo population at least out to ~ 85 kpc in this 
direction. On the other hand, it is difficult to discern whether an 
RGB sequence is present in the outermost bin (R > 65 kpc) along 
the minor axis. This asymmetry in the density of stars along the ma- 
jor and minor axes of the outermost elliptical bin suggests that the 
outer halo is either more flattened than the inner isophotes would 
suggest or that it is spatially inhomogeneous on the scales probed 
by VIMOS. The red extension of the RGB appears to become bluer 
with increasing radius indicating a decrease in the number of metal- 
rich stars. However, the colour separation between isochrones does 
not scale linearly with metallicity and we will subsequently see that 
this behaviour does not translate into a strong metallicity gradient 
(see Sect.O. 



4.2 Luminosity Functions as a Function of Radius 

In addition to CMDs, we also examine constraints on halo extent 
that come from consideration of the luminosity functions (LFs). We 
construct the LFs by selecting stars in a colour-magnitude box that 
minimizes the background/foreground contamination as much as 
possible. Stars are chosen in the magnitude range 23 < I < 25 
and in a diagonal colour range with the bottom edge comprising 
1 < (V - I) < 2 and top edge comprising 1.4 < (V - I) < 2.5 
(see Fig. [5]). We choose magnitude bins of ~ 0.3 mag, the largest 
photometric error within the adopted colour-magnitude range for 
the innermost radial bins, and derive histograms in the same radial 
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Figure 4. Top panel. Dereddened CMD for the combined four VIMOS 
pointings. Individual points are replaced by a Hess (density) diagram where 
the density is highest. We show the boxes used to select the RGB and AGB 
subsamples (red boxes). We overlay Dartmouth stellar isochrones for the 
RGB stars (green lines), with a fixed age (12 Gyr) and varying metallicities 
([Fe/H] = -2.5, -0.9 and -0.3). We also plot the p ost He-burning AGB 
phase of the Padova isochrones (cyan dashed lines, iGirardi et alj 120101) . 
The latter have ages of 2 and 4 Gyr and a metallicity of Z = 0.004 (see 
text for details). Central panel. CMD of all the extended objects rejected by 
quality cuts in our PSF photometry. A Hess diagram indicates the highest 
density region. Bottom panel. CMD of Galactic foreground stars as pre- 
dicted by the Besancon models (Robin et al. 2003, see text for details). The 
foreground simulations have been convolved with photometric errors and 
the observational incompleteness effects have been taken into account (see 
text for details). In the central and bottom panels we also overplot our RGB 
and AGB selection boxes to ease the comparison with the observed CMDs. 
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Figure 5. Dereddened Hess (density) diagrams as a function of elliptical radius, for both major (upper panels) and minor (lower panels) axes. We subdivide 
the stellar samples into three elliptical bins per axis (see also Fig. [I}. The colourbars (on a square root scale) indicate the number of stars per unit area and per 
0.05 x 0.05 mag 2 . We overlay Dartmouth stellar isochrones with a fixed age (12 Gyr) and varying metallicities ([Fe/H] = -2.5, -1.3, -0.7 and -0.3 dex). Note 
the different stellar densities among the radial bins. Finally, in the upper left panel we also show the selection box used to derive LFs (see Sect. 14.2b . 



bins as before. The LFs are corrected for incompleteness as a func- 
tion of magnitude and radius, and the expected foreground counts, 
as simulated by the Besancon models, are subtracted off. Finally, in 
order to facilitate comparison, the LFs in each bin are normalized 
to the total number of stars brighter than I = 24.6 (correspond- 
ing to the 50% completeness limit for the most crowded bin, i.e., 
R < 45 kpc along the major axis). 

Fig. [6] shows the LFs in radial bins along the major and minor 
axes. Poissonian errorbars are shown, which include the uncertainty 
in the completeness corrections and in the Besancon model counts. 
Although the LFs look broadly similar, the outer bins show an ex- 
cess of sources brighter than the TRGB (Iojrgb = 23.81 + 0.35) 
with respect to the innermost bin. Additionally, the LFs become 
progressively flatter at faint magnitudes in the outer radial bins 
with the variation being more pronounced along the minor axis. 
A Kolmogorov-Smirnov (KS) test rejects the null hypothesis that 
the LFs at radii beyond 45 kpc come from the same underlying 
population at a 99% significance level. 

Although these differences could result from genuine radial 



variations in the halo populations, it is also plausible that they re- 
sult from varying relative contributions of Cen A stars and contam- 
inants. Indeed, in the outer halo the latter will contribute an increas- 
ing fraction of the overall source counts and will thus have a more 
profound impact on the LF shape. In Fig. [7] we plot the LF per 
unit area of sources rejected as stars by our quality cuts, derived 
from the same selection box as above. If unresolved background 
contaminants follow the same luminosity distribution as resolved 
galaxies, then we can assume this is how these objects will con- 
tribute to the overall LF. We also show in Fig. [7] the LF per unit 
area of point sources in the 75 < R < 85 kpc radial bin along 
the minor axis, which has not been corrected for incompleteness 
or foreground subtracted (we use the bin 75 < R < 85 kpc rather 
than 65 < R < 85 kpc in order to avoid as much as possible the 
presence of Cen A halo stars). Despite the fact that the true frac- 
tion of contaminants remains unknown, this comparison illustrates 
the excellent agreement between the two estimates of background 
LF shapes at magnitudes brighter than the TRGB. This is reassur- 
ing as it indicates that the excess of bright stars with respect to the 
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Figure 6. Luminosity functions at different elliptical radii as indicated, for both major (left panel) and minor (right panel) axes. The chosen elliptical bins are 
the same as in Fig. [5] The LFs are derived from stars in a box with 23 < To < 25 and 1 < (V - I)o < 2 (bottom edge) and 1.4 < (V - I)o < 2.5 (top edge). They 
are normalized to the number of counts for Iq < 24.6 (vertical red dash-dotted lines, i.e., the 50% completeness limit for the R < 45 kpc subsample along the 
major axis). For each LF we report the number counts for Iq < 24.6. The LFs are corrected for incompleteness, and simulated foreground counts have been 
subtracted for each magnitude bin. The errorbars are poissonian. 
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Figure 7. Luminosity functions for the 75 < R < 85 kpc field along the 
minor axis (blue line) and for the extended sources rejected by our photo- 
metric cuts (black dashed line). The number counts are per unit area, and 
both curves are not corrected for incompleteness. 



innermost radial bins is very likely due to an increasing fraction 
of contaminants at large radii. At fainter magnitudes, the LF of re- 
solved sources turns over with respect to that of point sources in 
the 75 < R < 85 kpc minor axis bin, consistent with the progres- 
sive flattening seen in the outer LFs in Fig. [6] Clearly some of this 
behaviour will be a result of incompleteness in the resolved source 
counts but we have no way to correct for this. Without an accurate 
measurement of the shape and normalisation of the contaminant LF, 
it is impossible to quantify what fraction of the faint sources in the 
75 < R < 85 kpc bin LF are genuine Cen A halo stars however it is 
likely that this number is non-negligible. Indeed, we have derived 
TRGB values for each radial bin along both the major and minor 
axes and find all of them to be consistent with the value derived 
from the total stellar sample (see Sect.|3}, within the uncertainties; 
this supports the idea that a significant RGB population is present 
in all our fields. 



4.3 Radial Density Distribution 

In Fig. [8] we show the spatial distribution of sources falling within 
our RGB and AGB selection boxes. Small holes in the distribu- 
tion are mostly due to saturated stars. As noted before, the highest 
stellar density occurs in the maj 1 pointing. Additionally, Fig. [8] re- 
veals that the enhanced density in this region is due in part to a 
spatially-coherent substructure which is most prominent on the in- 
nermost chip; this overdensity is apparent in both the distribution 
of RGB and AGB candidates. At a radius of ~ 40 kpc, the new 
substructure lies well beyond the known shell system. To further 
explore this feature, we retrieved a scanned UK Schmidt Telescope 
(UKST) Ilia- J survey plate centred on Cen A and spanning 1.4° 
on a side. Fig. [9] shows a slightly smoothed version of this plate 
with an outline of the innermost chip of the majl pointing over- 
laid. The substructure detected in our star count analysis is clearly 
visible on this deep plate. Furthermore, the plate reveals that the 
substructure is part of a very large overdensity of stars which pro- 
trudes from the main body of Cen A along the north-eastern maj or 
axis. Given the morphology of this feature, the most likely origin 
is post-merger debris from a significant merger event and not indi- 
vidual dwarf galaxy accretions. 

Aside from the innermost major axis field, the RGB stellar 
density declines smoothly with radius along both axes, while for 
candidate AGB stars it remains relatively constant. We show the 
radial density profiles for the RGB sample in the left panel of Fig. 
ITOl These are computed using the same elliptical annuli as before. 
Star counts are calculated in radial bins for 30 < R < 90 kpc, in 
steps of 10 kpc, and corrected for incompleteness as a function of 
radius and magnitude. These profiles are not corrected for fore- 
ground/background contaminants. Instead, we indicate the RGB 
number density for the 75 < R < 85 kpc bin along the minor axis 
with a dash-dotted line, and assume that this represents an upper 
limit on the contamination. Similarly, a lower limit is given by the 
dashed line which is constructed from the sum of resolved back- 
ground source density and simulated foreground star density, de- 
rived in the same selection box as for the RGB sample. As can be 
seen, these contaminant estimates differ by a factor of ~ 4. Com- 
pleteness corrections are very significant for R < 40 kpc while the 
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Figure 8. Distribution on the sky (standard coordinates with respect to the center of Cen A) for the candidate RGB and AGB stars along the major (upper 
panels) and minor (lower panels) axes (see Fig.[T]for a reference). Red dashed ellipses are drawn in steps of 10 kpc, starting at an elliptical radius of 30 kpc out 
to 80 kpc. Black lines are, respectively, major and minor axes. Holes in the stellar distributions are mostly due to the presence of foreground saturated stars. 
Note that the highest density region is localized to the innermost chip along the major axis, for both RGB and AGB stars. To underline the highest density in 
the RGB sample, we replace individual points with a contour plot for densities larger than ~ 670 stars/arcmin 2 . 



last major axis point is calculated from only a very small area (see 
Fig©. It is clear that the surface density of RGB stars is always 
higher along the major axis than the minor axis, and especially for 
radii < 55 kpc. On the other hand, both axes exhibit similar profile 
shapes with steep declines in the number counts out to ~ 55 kpc 
and much flatter declines beyond that. 

Unfortunately, we cannot construct a surface density (or sur- 
face brightness) profile across the whole extent of Cen A's halo 
since the UKST plate does not have the required level of linearity 
for an accurate photometric calibration. Previous studies have con- 
strained the structure of Cen A over s mall a reas. Using data within 
the inner ~ 12 kpc, Ivan den Berghl dl976h fit a de Vaucouleurs' 



law ddeVaucouleurslll959h with parameters R ef [ = 6.1 kpc and 
V(#eff) = 22.15 mag arcsec" 2 . Additionally. iDufour et all dl979h 
find R eS = 5.5 kpc from U- and V-band profiles along the major 
axis, with V(R eS ) = 22.00 mag arcsec -2 . The smaller effective ra- 
dius found in the latter study probably results from their exclusion 
of data at radii > 4.7 kpc. 

It is of obvious interest to explore whether the stellar density 
profile in the outer halo behaves as expected from an extrapolation 
of the inner profile and we attempt to address this as follows. In the 
right panel of Fig. [TO] we present contaminant- subtracted versions 
of our RGB radial density profiles using both the upper (filled sym- 
bols) and lower (open symbols) limit estimates for the contaminant 
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Figure 10. Left panel. Radial density profiles for RGB stars along both major (filled black circles) and minor (filled light blue triangles) axes, as a function 
of elliptical radius. The radius is also expressed in terms of effective radii (top axis). Star counts have been corrected for incompleteness (both radial and in 
magnitude). The dashed-dotted red line indicates the upper limit estimate for background/foreground contamination (i.e., the RGB number density for the 
75 < R < 85 kpc bin along the minor axis), while the dashed blue line represents the lower limit estimate (i.e., the number density of resolved background 
sources and simulated foreground stars). Right panel. Same as the left panel, after subtraction of the upper limit estimate (filled symbols) and lower limit 
estimate (open symbols) for background/foreground. Overplotted for each axis and for each of the two subtracted profiles (solid lines for the upper limit and 
dashed lines for the lower limit for background/foreground) is a de Vaucouleurs' profile with effective radius of 6.1 kpc and arbitrarily scaled to match the 
datapoint at R = 55 kpc. The right-hand vertical axis reports surface brightness values corresponding to the de Vaucouleurs' profile arbitrarily scaled to the 
(upper-limit) contaminant- subtracted major axis profile. The errorbars are poissonian. 




Figure 9. Extended UKST photographic plate image (see also Fig.[T] north 
is up, east is left). Overlaid is the location of the innermost VIMOS chip 
along the major axis (red rectangle). 

le vel. We overplot the de -reddened de Vaucouleurs' profile derived 
bv lvan den Berghl(ll976h . arbitrarily scaled to match the R = 55 kpc 
datapoint in each case. Our VIMOS density profiles appear broadly 
consistent with fall-off of the extrapolated de Vaucouleurs' profile. 
There is a reasonable agreement along the minor axis when the 
upper limit on the contaminant level is adopted, while using the 
lower limit yields an excess of stars at radii R > 65 kpc. The ma- 
jor axis profile deviates significantly from the inner extrapolation 
at R > 75 kpc regardless of what contaminant level is subtracted. 
Although some part of this behaviour could be due to a genuine flat- 



tening of the halo profile at these radii, the increased radial density 
in the R ~ 85 kpc bin along the major axis suggests additional low- 
level substructure in these parts. Visual inspection of the outermost 
chip along the major axis (Fig. [8} provides a tantalising suggestion 
of this, with a marginally higher stellar density in the outer portion 
of the chip. Had we assumed a de Vauco uleurs' profile with a n ef- 
fective radius of 5.5 kpc, as derived by Dufouretal our 
conclusions would be unchanged. For reference, in the right panel 
of Fig. \\0\ we indicate the surface brightness values correspond- 
ing to the extrapolated de Vaucouleurs' profile, arbitrarily scaled to 
match the (upper-limit) contaminant- subtracted major axis profile 
at 55 kpc. 

Finally, the radial density profiles for candidate AGB stars (not 
shown here) are essentially flat along both axes, except for the sig- 
nificant major axis overdensity at radii < 55 kpc. Any genuine AGB 
stars in these parts are most likely completely outnumbered by fore- 
ground and background contaminants. 



5 METALLICITY DISTRIBUTION FUNCTIONS 

For largely coeval populations, the colour and width of the RGB 
constrains the mean metallicity and metallicity spread of the popu- 
lation (the RGB colour is m ore sensitive to cha nges in m etallicity 
than it is to changes in age. IVandenBerg et al.l e.g.. 120061 .). Given 
that we have not found evidence for a significant intermediate- age 
AGB population in our VIMOS fields (except for the innermost 
radial bin along the major axis), we will assume that the Cen A 
halo population is essenti ally coeval and p r oceed to derive photo- 
metric MDFs (see e.g., lHarris et al.|[l999l: lHarris & Harrisll2002l : 
iReikuba et al J20Q5L and references therein). 

We again adopt 12 Gyr or-enh anced ([cr/Fe] = +0.2) 
isochrones from the Dartmouth group dDotter et al.l 12008k with 
metallicity ranging from [Fe/H] = -2.5 to [Fe/H] = -0.3, sam- 
pled every 0.2 dex. This choice of age is driven by the analysis of 
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Figure 11. Metallicity distribution functions as a function of elliptical radius, for both major (left panels) and minor (right panels) axes (solid black lines). 
The MDFs are derived via isochrone interpolation with a fixed age (12 Gyr) and varying metallicity. The histograms are corrected for incompleteness and 
normalized per unit area, and the median metallicity (after completeness correction) is reported for each panel. The errorbars are computed via MC simulations, 
taking into account photometric errors and random errors stemming from the interpolation technique. Finally, in each subplot we also show the MDF for the 
75 < R < 85 kpc bin along the minor axis (red dot-dashed histograms, i.e., upper limit estimate for the foreground/background contamination), and the 
combined MDF for resolved background sources and simulated foreground stars (blue dashed histograms, i.e., lower limit estimate). 



Reik uba et al.1 (1201 lh which indicates that the bulk of the star for- 
mation in their 38 kpc field occured ~ 12 Gyr ago on a relatively 
short timescale (~ 3 Gyr). 

Our MDFs are calculated only using stars which fall within 
our RGB selection box, to ensure good completeness. In this region 
of the CMD, the stellar isochrones are also more widely separated 
thus reducing the uncertainty in individual metallicity determina- 
tions. We linearly interpolate within the isochrone grid to deter- 
mine the metallicity of each star. The uncertainties are computed 
by producing 1000 Monte-Carlo (MC) realizations of the interpo- 
lation process: each time the CMD position of the star is varied 
within its photometric uncertainty (with a Gaussian distribution), 
and the metallicity recomputed. 

While the derivation of precise absolute metallicity estimates 
is not possible with photometric data alone, we can nevertheless 
get robust results on the relative metallicities within our sample. 
Fig. [TTI shows the MDFs, corrected for incompleteness, as a func- 
tion of radius along both axes, while Tab. [3] tabulates their main 
properties. The MDFs exhibit the same overall shape consisting 
of a peak around [Fe/H]~ -0.8 and a long tail to lower metallic- 



ities. At R < 65 kpc, the major axis has a slightly larger num- 
ber of stars more metal-rich than [Fe/H]~ -0.70 compared to the 
minor axis, which translates into a small increase in the median 
metallicity. This mild enhancement could plausibly be related to 
the presence of the substructure we have uncovered in these parts. 
Intriguingly, we find the median [Fe/H] is rather constant across 
the radial extent of our survey, varying from only -0.88 to -1.04. 
In particular, there is a change of A [Fe/H] = -0.16 dex along the 
major axis (~ 35 - 80 kpc, see Fig [8} corresponding to a gra- 
dient of ~ -0.004 dex/kpc. Along the minor axis, the change is 
A [Fe/H] = -0.1 dex over ~ 30 - 85 kpc and thus a similarly small 
gradient of ~ -0.002 dex/kpc. The metallicity spread, as estimated 
by the standard deviation of the individual metallicity measure- 
ments, is also essentially constant at ~ 0.45 dex throughout the 
halo. Finally, the fraction of metal-poor stars (which we define 
as f mp = N[F e /#]<-i.o/Nfof) shows only a very modest increase of 
~ 10% with radius, reaching a maximum of ~ 50% in the outer- 
most bin. 
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Table 3. Median metallicities and metal-poor stars fractions (f mp = ^[Fe/H]<-i.of^tot) as a function of radius, for different background/foreground 
subtractions. 



Field 


[Fe/H] med 


\^&[^\med,sub-upp^ 


[pQ/H\med,sub-low^ 


fmp 


fmp,sub-upp 


fmp, sub-low 


R < 45 kpc (major) 


-0.88±0.04 


-0.87±0.04 


-0.87±0.04 


0.40±0.04 


0.39±0.05 


0.39±0.05 


45 kpc < R < 65 kpc (major) 


-0.89±0.01 


-0.85±0.01 


-0.87±0.01 


0.41±0.01 


0.34±0.01 


0.38±0.01 
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a after subtraction of the upper limit on background/foreground estimate, i.e., the MDF of the 75 < R < 85 kpc bin along the minor axis. 
^ after subtraction of the lower limit estimate, i.e., the combined MDF of resolved background sources and simulated foreground stars. 
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Figure 12. Normalized metallicity distribution functions for the combined 
VIMOS fields at 35 < R < 45 kpc along major and minor axes (black solid 
line), compared to the one for the outermost ACS/HST pointing (red solid 
line). The solid line MDFs are drawn from stars in the RGB selection box 
shown in Fig. [4] and the VIMOS MDFs have been corrected for incomplete- 
ness. The red dot-dashed MDF is derived using all RGB stars brighter than 
Iq = 26.3 and within the whole colour range of the ACS/HST photometry. 
The VIMOS errorbars are drawn from MC simulations (see text) , while the 
ACS errorbars are poissonian. 



5.1 Uncertainties 

There are a number of uncertainties in the derivation of the MDFs 
that we now consider. In particular, we wish to address to what ex- 
tent these uncertainties could affect our derived MDFs as a function 
of radius. Firstly, intermediate- age AGB stars ascending the giant 
branch could contaminate our RGB samples and this could be im- 
portant in our innermost major axis chip where the clear signature 
of such a population exists. AGB stars will lie on the blue side of 
the RGB, and will thus artificially enhance the metal-poor popula- 
tion in our MDFs. However, given that the metal-poor population 
in the innermost major axis bin is small (see Fig.[TTJ, we conclude 
that this cannot be a significant effect. 

Secondly, similar to the LFs discussed in Sec. 14.21 the MDFs 
we derive will contain both genuine Cen A halo stars as well as con- 
taminants. As the contaminants are not RGB stars at the distance of 
Cen A, the metallicities computed for them will be meaningless. 
Both the number and the colour distribution of the contaminants 



can potentially influence our results. In Fig. Qj] we show normal- 
ized "MDFs" for our two estimates of the contaminant population, 
derived in exactly the same way as described above. Our lower 
limit estimate, based on summing resolved galaxies and foreground 
models, yields an "MDF" that is fairly flat while our upper limit es- 
timate, based on the outermost portion of the outer minor axis field, 
bears a strong similarity to those we have derived for the halo. This 
is not surprising since this field may still contain a significant num- 
ber of genuine Cen A RGB stars. We subtract these estimates of the 
contaminant "MDF" from our measured MDF and recompute the 
median metallicity and metal-poor fractions (see Tab. [3}. This has 
only a small impact on our results and does not affect our inferences 
regarding the metallicity gradient and metal-poor fraction. 

Finally, although we have corrected for incompleteness in de- 
riving our MDFs, these corrections are very significant for the red- 
dest ((V - I)o > 2.5) populations that fall within our RGB selection 
box. In order test how accurate our completeness corrections are 
for these metal-rich populations, we use the photometric catalogue 
for the HST/ACS fie ld at 38 kpc, kindly provided by M. Rejkuba 
(Reik uba et al.| [2005), which reaches to redder colours and fainter 
magnitudes than our VIMOS data. We deredden the photometry 
with the mean extinction values reported in Sect. For the magni- 
tude and colour ranges of relevance for comparing to our VIMOS 
data, the photometric and colour uncertainties are generally smaller 
than ~ 0.03 and ~ 0.04 mag respectively, and the completeness as 
derived from artificial star tests is close to ~ 100%. In addition to 
photometric accuracy, the other main advantage of this catalogue 
over our ground-based data is its purity. The high resolution (0.049 
arcsec/pixel) of the ACS means that background galaxy contami- 
nants can easily be identified and rejected. 

We initially derive the MDF for the HST/ACS catalogue us- 
ing the same RGB selection criteria and procedure as for our VI- 
MOS data. In Fig. Q21 (solid lines), we show the resulting ACS 
MDF as well as the VIMOS MDF for the same radial range 
(35 < R < 45 kpc), constructed using data along both major and 
minor axes. The HST/ACS sample yields a median metallicity of 
[Fe/H] = -0.89, and a spread (standard deviation) of ~ 0.44 dex. 
For a comparison, the VIMOS MDFs at these radii have a median 
metallicity of [Fe/H] = -0.90 with a spread of ~ 0.42 dex. Reassur- 
ingly, there are only small differences in the median metallicity and 
the shape of the two MDFs, confirming that that our completeness 
corrections for the reddest stars in the VIMOS data are accurate, 
and that contaminating galaxies have a negligible influence on our 
results at these radii. 

We then proceed to derive the MDF for the HST/ACS data 
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using a broader range of colour and magnitude than we are sen- 
sitive to with VIMOS in order to assess what systematic biases 
might arise from the limitations of our photometry. We select RGB 
stars brighter than I = 26.3 and use the whole range of colours 
in the HST/ACS CMD. The dashed line in Fig. OH shows the re- 
sultant MDF in this case. Overall, there is a good agreement be- 
tween this MDF and our VIMOS MDF however it peaks at some- 
what higher metallicities and has a metal-rich tail not seen in the 
ground-based data. The median metallicity derived in this way is 
[Fe/H] = -0.75 with a spread of 0.46 dex. This value translates into 
[M/H] = -0.61 for our choice of or-enhancement, and is in excellent 
agreement with the mean metallicity [M/ H] = -0.65 reported by 
iReikuba et all 12005) (obtained using the IVandenBerg et al.ll2000l 
evolutionary models). Finally, the metal-poor fractions for the 
HST/ACS MDFs are 0.38 (VIMOS selection box) and 0.29 (full 
colour range) which are within ~ 10% of our VIMOS estimate 
(i.e., 0.40 for the considered radial range). 

In summary, none of the effects discussed above are expected 
to have a significant impact on our results. The lack of sensitivity to 
very red stars in the VIMOS data may slightly bias (~ 0. 15 dex) our 
metallicity estimates towards lower values, but it is unlikely that 
this effect could disguise any strong radial gradients. Contamina- 
tion from background/foreground objects has even less of an effect 
on our measurements. We conclude that the moderately high and 
relatively uniform metallicity in the extended halo of Cen A is a 
robust result. 



6 DISCUSSION 

6.1 Halo Extent and Structure 

The intrinsic faintness of the outskirts of gEs has made it challeng- 
ing to quantify the structure and extent of these objects, even in our 
local universe. We have shown the existence of an RGB popula- 
tion to at least ~ 85 kpc (~ 14R eS ) along the major axis of Cen A. 
Along the minor axis, there is a clear RGB signature in the CMDs 
to ~ 65 kpc, however analysis of LFs and TRGB distances sug- 
gests that the stellar population extends further in this direction too. 
The RGB population in Cen A is therefore at least as extended as 
the distribution of planetary nebulae and GCs JPeng et al.ll2004al lbl; 
IWoodlev et al. 2007, 2010b; Har ris et alJ2012h . A few studies have 
show n the existence of stars at very large r adi i in other ellipticals 
(e.g. IWeil et all 1 1991 lHarris et all 12007 iTal & van Dokkuml 
l201ll) . Taken together, these studies demonstrate the vast extents 
of elliptical galaxies. Almost all previous studies of ellipticals have 
focused on regions lying at < 1 - 2R eS . While this is sufficient 
to sample the bulk of the light in these systems, it samples only 
a small fraction of their physical extents. With the increasing fo- 
cus o n questions related to the size evolution of early-type galaxies 
(e.g., iBuitrago et all boOS), an understanding of the properties and 
nature of the peripheral regions of ellipticals is of paramount im- 
portance. 

The RGB stellar density in Cen As halo is always higher along 
the major axis than the minor axis for a given elliptical annulus. 
Within ~ 55 kpc (< 9R eS ), this results from the presence of a co- 
herent overdensity of stars that we have uncovered along the major 
axis. The morphology of this feature suggests it is post-merger de- 
bris, perhaps originating from the same recent (< 1 Gyr) ev ent tha t 
led to the galaxy's peculiar inner structure (see, e.g., Ilsraelll 19981) . 
The feature is unlikely to be due to triggered star formation from 
Cen As radio jets given the estimated timescale (10 7-8 years) for 



this process compared to the typical > 1 Gyr age of RGB stars 
dlsraellll998h . The existence of such substructure in Cen As far 
outer halo indicates that these regions are significantly inhomoge- 
neous and that a wide-field approach is required to properly study 
and interpret the resolved stellar populations in these parts. 

Beyond 55 kpc, the difference between the major and minor 
axis density profiles could plausibly be explained if our assump- 
tion of elliptical isophotes with b/a=0.77 and a position angle of 
35° breaks down. Indeed, the two profiles could be brought into 
better agreement if b/a~0.5-0.6 in these parts, implying a higher el- 
lipticity in the outer halo. Interestingly, ITal & van Do kkum (2011) 
find the same result in their stacking analysis of SDSS luminous red 
galaxies. In particular, they find a mild increase in ellipticity from 
0.25 to 0.3 (or equivalently b/a varying from 0.75 to 0.7) at radii 
> 5Rqs . On the other hand, the profile differences in Cen A could 
simply be due to additional low-level substructure further out along 
the major axis, a tantalising hint of which may exist in our data. 

Finally, we have examined whether the star count measure- 
ments we have made in the outer halo are consistent with an ex- 
trapolation of the inner surface brightness profile. Along both the 
major and minor axes, there is evidence for a flattening in the pro- 
file shape beyond ~60-70 kpc (> 10-11 .5R eS ), the exact details of 
which depend on the contaminant level subtracted and how the pro- 
files are normalised. Subtracting the upper limit on the contaminant 
level, approximate power-law profiles of R~ 3 3 and R~ 3 8 can be fit to 
the major and minor axis profiles of Cen A, respectively. A flatten- 
ing of the outer profile was also seen beyond ~ 8ff e ff in the stacked 
surface brightness profile of ITal & van Dokkuml (l201lh . These au- 
thors were not able to tell whether the excess light at large radius 
was due to a flattening in the galaxy surface brightness profile, or 
to intragroup light. Indeed, the same ambiguity persists to some ex- 
tent here although the Cen A group is unlikely to have an extensive 
intragroup light component. 

6.2 Metallicity Gradient 

The photometric MDFs derived from our VIMOS data yield me- 
dian values of [Fe/H] me j ~ -0.9 to -1.0 dex which, given our 
choice of ^-enhancement, correspond to [MfH] med ~ -0.75 to 
-0.85 dex. The median metallicity is essentially constant, varying 
by only ~ 0.1 - 0.15 dex across the radial extent of our survey 
(~ 5 through 14R e s). This variation implies negligible gradients, 
< -0.004 dex/kpc. There are also no strong azimuthal variations, 
aside from a small enhancement in metallicity in the region of the 
outer halo major axis substructure. This suggests that either the stel- 
lar populations of the accreted system do not differ substantially 
from the native populations in the outer halo of Cen A or that sig- 
nificant mixing of the debris has already occurred. 

The moderately high metallicity found argues against a sce- 
nario in which most of the outer halo stars have been accreted 
from very low mass objects. ICrnoievic et al.l (l2010h study a sam- 
ple of dwarf satellites of Cen A with luminosities in the range 
M v ~ -10.7 to -13.9, corresponding to ~ 0.005 - 0.1% the lumi- 
nosity of Cen A. They find median metallicities of [M/H] < -1.0, 
considerably below the outer halo metallicity of Cen A How ever, 
this result may not be surprising. Indeed, lOser et all d2012l) use 
simulations to argue that the dominant mode of halo assembly is 
through minor mergers with typical mass ratios of « 1:5; such sys- 
tems are likely to be more metal-enriched than low mass dwarfs. 

We can compare our results to previous HST studies of 
Cen As halo. These cover 4 fields in the radial range fro m ~ 1.3 to 
~ 6.2/^eff, none of which lie along the major axis (Harri s et al.ll999l: 
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lHarris & Harrisll2000l 120021: iReikuba et aDl2005h . Aside from the 
central pointing which contains a combination of halo and bulge 
stars, these fields all exhibit a moderately high mean metallicity 
([M/H]~ -0.4 to -0.6) with only a small population (~ 10%) of 
low metallicity stars. When combined with our results presented 
here, this suggests a variation of only < 0.5 dex in [M/H] over 
the entire extent of the halo. Harris et al. were able to reproduce 
their MDFs with a simple chemical evolution model in which early 
star formation proceeds with an initial stage of rapid infall of very 
metal- poor gas, after which t he infall dies away exponentially. Sim- 
ilarity, Rejkuba et al. bonh argue that most of the stars in their 
field at ~ 6.2R eff formed rapidly and at very early times, with a mi- 
nority intermediate- age component having formed more recently. 
The moderate metallicities we find at larger radii in the halo sug- 
gests that the qualitative conclusion of rapid collapse continues to 
hold in these parts, although we have yet to model our MDFs in de- 
tail. Interestingly, recent kinematic and abundance studies of GCs 
in Cen A's halo (out to ~ 40 kpc) also support a fa st, early forma- 
tion ph ase, coupled with subsequent merger events dWoodlev et al.l 
l2010al lbh. 

Unfortunately, it is not straightforward to compare our results 
to those for other gEs. Very few observational studies have been 
able to investigate the metallicity content of stellar populations at 
galactocentric radii as large as we have probed in Cen A. HST/ACS 
observations have been obtained centre d at ~ l2R eS for the nearby 
"classical" gE NGC3379 (~ 11 Mpc: lHarris et al.ll2007ah and at 
~ 3/? e ff f° r the inter mediate- sized elliptical, NGC3377 (~ 10 Mpc; 
Harris et al . 2007b). The resulting picture is different for these 
two systems. NGC3379 shows an extremely broad and flat MDF, 
with the metal-poor population becoming increasingly dominant 
at > lOfleff • On the other hand, NGC3377 has an MDF peaked at 
~ -0.6 dex and does not show a significant gradient out to the 
furthest radius probed (~ 5R eS ). In integrated light analyses, some 
systems show cont i nuing metallicity gra dients out to ~ 2 - 4R eii 
dGreene et al .120121 ; I Weii mans et al .12 009) while o thers show steep- 
ening of gradients when traced a s far out as 8/? e ff (L a Barbera et al.l 
120121) . The stacking analysis of iTal & van Dokkuml (1201 ll) yields 
a clear colour gradient within the inner ~ 3i? eff , with the popula- 
tions getting bluer and thus more metal-poor (and possibly older) 
with radius, while the color index flattens out at larger galactocen- 
tric distances (out to ~ 14i? e ff)- Larger samples will be required to 
establish what is the typical behaviour for the metallicity gradient 
at large radii in gEs. 

Comparison to the predictions of theoretical models is sim- 
ilarly difficult since existing models do not explicitly address the 
large radii probed by our study. Theoretical models of gE formation 
are able to predict both the presence and the absence of metallic- 
ity gradients at smaller radii. Monolithic collapse models, in which 
gas sinks to the centre of the potential well and is enriched by the 
first generations of evolving stars, predict steep radial metallicity 
gradients which become more pron ounced with inc reasing g alaxy 
mass due to gas loss efficiency (e.g.jMatteucciiri984l ; IChiosill2002l ; 
iKawata & Gibsoijl2003l ; lKobavashi|[2004l ). On the contrary, in hi- 
erarchical models, the physical properties of the resultant galaxy 
will largely depend on those of the progenitors however the expec- 
tation is for metallicity gradients to be shall ower due to the mixing 
of stars (Kobayashi 2004; Naab et al. 2009). Whether or not these 
inner trends should extend to larger radii remains an open question. 

Finally, we note in passing that the metallicity of RGB stars 
in Cen A's outer halo is very similar to that which we have 
foun d in the extended stellar envelopes around some spiral galax- 
ies (Barker et al. 2009, 2012). These extended envelopes are also 



characterised by flatter surface brightness profiles than seen in the 
inner regions. It is interesting to speculate whether moderate metal- 
licity extended stellar structures could be an ubiquitous feature of 
all large galaxies, regardless of morphological type. 



7 CONCLUSIONS 

We have conducted a wide-field survey of resolved stellar popula- 
tions in the outer halo of our closest gE, Cen A, using the VIMOS 
imager mounted on the VLT. Two fields were imaged along each 
of the major and minor axes sampling projected elliptical radii in 
the range - 30 - 85 kpc (or ~ 5 - 14tf eff ). From PSF-fitting pho- 
tometry, we derive CMDs in V- and /-bands that reach ~ 2 mag 
below the TRGB. The CMDs show evidence for an old RGB pop- 
ulation, which can be traced to the furthest extent of our survey 
(~ 85 kpc) along the major axis and to at least ~ 65 kpc along 
the minor axis. However, consideration of LFs and MDFs strongly 
suggests that the stellar population along the minor axis extends to 
at least ~ 85 kpc too. 

The spatial distribution of RGB stars has been mapped as a 
function of radius and azimuth. We have uncovered a prominent 
localized overdensity extending to ~ 55 kpc (~ 9R e s) along the ma- 
jor axis. Inspection of a deep UKST plate confirms the existence of 
this structure and shows a morphology consistent with post-merger 
debris. Our VIMOS data indicate that the metallicity of the con- 
stituent stars in this structure do not differ significantly from that 
of the underlying halo population. Even beyond the extent of this 
substructure, the major axis stellar density is higher than the minor 
axis one at a given elliptical radius. This could be evidence for an 
increasing ellipticity in the outer halo, or conversely for residual 
low level debris contaminating the outer halo major axis. Uncer- 
tainties in the contaminant level make it difficult to rigorously as- 
sess whether the outer halo density profile follows an extrapolation 
of de Vaucouleurs law which characterises the inner regions. Both 
the major and minor axis profiles appear generally consistent with 
such an extrapolation, but there is evidence for a flattening in the 
profile beyond - 60-70 kpc (~ 10 - 11.5fl e ff). 

We derive photometric MDFs via isochrone interpolation of 
RGB stars. The median metallicities we find are relatively high 
(<[Fe/H]> meJ ~ -0.9 to -1.0 dex, corresponding to [M/H] ~ -0.75 
to -0.85 dex), with broad spreads of ~ 0.45 dex. We analyse the 
median metallicities as a function of radius in order to search for 
radial gradients, but find only a modest decrease of ~ 0.1 -0.15 dex 
over a range of ~ 5 to !4R eff . Moreover, the fraction of stars more 
metal-poor than [Fe/H] = -1.0 increases by only ~ 10% as a func- 
tion of radius, reaching ~ 50% at the outermost radii probed. When 
combined with previously published results, the stellar metallicity 
in Cen A can be seen to vary by < 0.5 dex over the entire radial 
range of 85 kpc. This apparent constancy of the metallicity is vey 
intriguing, and suggests that if Cen A does have a metal-poor halo 
component then it can only dominate much further out. 

Our results demonstrate that important constraints on elliptical 
galaxy formation can come from studying the extent and properties 
of their extended halos. The faintness and inhomogeneity of these 
parts requires an approach that combines both depth and wide -field 
coverage. Unfortunately, Cen A is the only gE that can, at present, 
be readily resolved into stars over large areas from the ground; 
other techniques, perhaps less optimal, will be required for addi- 
tional gEs. Our findings indicate that Cen A itself is also worthy of 
further study; observations at even larger radii than we have probed 
here will be required in order to establish its true extent and search 
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for evidence of an underlying metal-poor component to its stellar 
halo. 
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